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We derive the invariant measure on the manifold of multimode quantum Gaussian states, induced by the 
Haar measure on the group of Gaussian unitary transformations. To this end, by introducing a bipartition of 
the system in two disjoint subsystems, we use a parameterization highlighting the role of nonlocal degrees of 
' freedom - the symplectic eigenvalues - which characterize quantum entanglement across the given bipartition. 

A finite measure is then obtained by imposing a physically motivated energy constraint. By averaging over the 
• local degrees of freedom we finally derive the invariant distribution of the symplectic eigenvalues in some cases 

^SJ ■ of particular interest for applications in quantum optics and quantum information. 



I. INTRODUCTION 



The development of quantum information theory has stimulated and motivated the study of the un- 
derlying geometry of quantum states and operations [6]. In particular, the analysis of the geometry 
of states describing composite quantum systems is fundamental because of the consequences of local 
and nonlocal transformations for the properties of correlations and entanglement [4, 5, 15]. A typical 
problem in quantum information is the characterization of the entanglement features across a given 
. partition of a quantum system into two (or more) subsystems. As a matter of fact, entanglement is the 

■ most important resource for quantum computation and quantum information. Nevertheless, it is not 

^ , yet completely understood. Therefore, it is of utmost importance to obtain normal forms for bipartite 

O^- quantum states, in which the role of the nonlocal degrees of freedom, characterizing the entanglement, 

is highlighted. For the case of a quantum system described by a Hilbert space of finite dimension 
n ~ C", and a bipartition 'H='Ha®'Hb with Ua ^ €"^^,^.3 ^ and + ris = n, a normal 
^ ' form of fundamental importance is given by the Schmidt decomposition [26] of a normalized vector 

. G H, representing a pure state of the system: 

in 

<N: IV') = 2Zv^li)A(^|j)B 

o: 

, Here {\j)A}j=i,...,nA^ {\j) B}j=i....,nB orthonormal systems of the local Hilbert spaces Ha and 

Hb, respectively, and we have assumed, with no loss of generality, that ua < ns. All states can be 
reduced to such a form by acting with local unitary transformations on the subsystems. The components 
of the probability vector p = {pi, . . . ,pn,^) E A„^_i, where A„^_i is the (n^ — l)-simplex (p^ > 0, 
' J2Pk ~ 1)' ^6 known as the Schmidt coefficients. They completely characterize the nonlocal features 

^ . of entanglement in pure states [35]. 

In the study of entanglement in bipartite systems, also in view of applications [16], one is interested 
in the typical features of entanglement [17], where the notion of typicality is defined according to a 
suitable probability measure. In particular, when dealing with finite-dimensional Hilbert spaces, one 
is interested in the distribution of the Schmidt coefficients. Although a certain arbitrariness exists in 
the choice of this distribution, the natural and most unbiased distribution is the one induced by the 
Haar measure on the relevant unitary group. For a system with Hilbert space H ~ C", this is the 
71-dimensional unitary group U(7i). The corresponding induced probability measure of the Schmidt 
coefficients is [19, 20] 

riA nA 

dfisip)^Pip)dp=a,,nA n iPh-PkfY[p]^-''^dp, (1) 

h>k=l j=l 
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with dp the Lebesgue measure on A„^_i, and Cn,nA ^ normaHzation factor. The probabiUty measure 
P{p) is by construction invariant under the action of the unitary group U(n), that is, it is independent 
of the reference frame. This expression has been the starting point of several studies concerning the 
typicality of entanglement, see e.g. [10-12, 14, 21, 24]. We stress that the derivation of a measure for 
a (physically relevant) class of states is fundamental for the study and comprehension of the features 
of entanglement. Moreover, these kinds of studies can shed light on interesting details of the geometric 
structure of the Hilbert spaces used to describe quantum systems and the properties of unitary groups 
describing quantum dynamics. 

In the case of Hilbert spaces with infinite dimensions H = L^(R)®" ~ i^(R"), the limit of the 
invariant distribution (1) is not well defined. Then, in order to obtain a well defined distribution, and 
mostly for their relevance in physical applications, here we consider a specific finite-dimensional, yet 
unbounded, manifold Q C H of pure states, known as Gaussian states [9, 32-34, 36]. 

Entanglement typicality in multimode Gaussian states has been the subject of previous investiga- 
tions [1, 29-31], whose starting points were the probability measures induced by microcanonical-like 
and canonical-like ensembles and the concept of "concentration of measures". Here we follow a dif- 
ferent route and derive a different probability measure. More specifically, by merely relying on the 
symmetries of the manifold of Gaussian states, we derive the invariant measure induced by the Haar 
measure on the group of Gaussian unitary transformations, whose orbit passing through the vacuum 
state is the manifold of Gaussian states, see e.g. [33], and [22]. 

Furthermore, in view of physical applications, we introduce a suitable effective cutoff by means of an 
energy constraint that enables us to normalize the invariant measure, then to compute the average over 
the local degrees of freedom and finally to derive the probability distribution of the nonlocal parameters 
(symplectic eigenvalues) characterizing bipartite entanglement in Gaussian states. As an example, we 
consider a specific submanifold of Gaussian states which is most relevant for applications in the domain 
of quantum optics. 

II. DEFINITIONS AND MAIN RESULTS 

A. Multimode Gaussian States 

In this Subsection we briefly recall few basic definitions and properties of Gaussian states. For a 
comprehensive review of their mathematical features and physical relevance we refer to [9, 32-34, 36]. 

Definition 1 A continuous variable (CV) quantum system is defined by a collection ofn bosonic modes 
with canonical annihilation and creation operators 

ak,al k = l,...,n, 
obeying the canonical commutation relations 

[ah,a\]=5hk, [ah,ak\ = [a\,a\]=Q V/i, fc = l...n, 

where ^ denotes the adjoint operator The canonical operators act on the Hilbert space H = 
span{|{j'i, j2, . . . , j„})}js.=o,i....,oo> where the vacuum state |0) := |{0, 0, . . . , 0}) is characterized 
by afc|0) = for all k = 1, . . . ,n, and 

k=l 

In the following we will study the action of multimode Gaussian unitaries, namely those unitary oper- 
ators generated by quadratic self-adjoint polynomials of the canonical operators. 

Definition 2 Let us consider a Hamiltonian, quadratic in the canonical operators. 



Hg — Hf) + Hi + H2, 



(2) 



3 



with 

Ho = 0, (3) 

n 

Hi = iE(^fc4-e:«^) : (4) 



k=l 



H2 = Yl (Mhkahau + MllkO-Wk + ^hkalak) , (5) 



h,k=l 

where G M, ^ S C" , Af is a complex and symmetric matrix, N is a Hermitian matrix, and i denotes the 
imaginary unit. We define as multimode Gaussian unitary for n-mode Gaussian unitary) the exponential 
ofHa: 

UG = e-'"^. (6) 

The set of n-mode Gaussian unitaries is a group, denoted as the group of n-mode Gaussian unitaries. 
The quadratic Hamihonians of the form (2) constitute the Lie algebra of this group. It is well known that 
the group of n-mode Gaussian unitaries is a projective representation of the inhomogeneous symplectic 
group ISp(2?i, M), with the quadratic terms (5) corresponding to the homogeneous subgroup Sp(2?i, M), 
and the linear terms (4) to the subgroup of translations T2n, isomorphic to C". We hence denote as 
homogeneous Gaussian unitaries the Gaussian unitary transformations obtained by putting ^ = 0. 

Definition 3 We call any vector state of the form 

IV'g) =Z^g|0) =c"'^^|0), 

for some multimode Gaussian unitary IAg, a Gaussian pure state. The set of such states will be denoted 
hyQ. 

The manifold of Gaussian pure states C "H is hence the orbit of the action of the group of Gaussian 
unitary operators on the vacuum state. The orbit of the subgroup of homogeneous Gaussian unitaries is 
a submanifold of Gaussian states, that we refer to as the submanifold of homogeneous Gaussian states. 
Homogeneous Gaussian states are characterized by the conditions 

(V'GlafclV'G) = (V'GlallV'G) = 0, fc = l,...,n. 

According to the definition of Perelomov [25] Gaussian states are coherent states for ISp(27i, M), and 
the homogenous Gaussian states are coherent states for Sp(2n, R). 

Let us hence consider a bipartition of the system into two disjoint subsets, denoted A and B, of the 
canonical modes, 

aAkiO-jik , fc = 1, . . . , 
aBk,aBl , k = 1, . . . ,nB , 

with UA + ns = n. This naturally induces on the total Hilbert space the bipartition into two disjoint 
local Hilbert spaces, H = Ha <E) Hb- Without loss of generality we assume ua < tib. For Gaussian 
pure states of a bipartite system the following Proposition holds: 

Proposition 1 ([7]) By acting with local Gaussian unitary transformations 14 q andUQ on subsystems 
A and B, any Gaussian pure state |'0g) of ua + tlb modes can be put into the canonical form 

riA 

= cxp (^rfefl^^-as/c - rfeaAl.aB|,^|0) , (7) 
fc=i 

where ri, r2, . . . , are real and non-negative parameters {unique up to permutations) associated to 
the state IV'g)- 
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In order to clarify our notation, let us consider the case n = 2, with ua = = 1- Then any two 
mode Gaussian state, up to local Gaussian unitaries, can be put in the form 

|TMSV) = exp {raAttB - raA^as^) |0) , (8) 

which is the so-called two-mode squeezed vacuum (or twin-beam state), a state of utmost importance 
in quantum optics and in CV implementations of quantum information processing [9]. 
By expanding the exponential in Eq. (8) one gets 

I — 9 — °° /u -^\ 

|TMSV)^y-^g(-^) 

where 

|j) = (j!)-V2(«t).|o), 

and V = cosh 2r. By taking the partial trace over one of the two modes, one gets a thermal-like reduced 
state 




The entanglement of the two-mode squeezed vacuum is hence characterized by the single parameter 

ly e [l,+oo). 

Coming back to the multimode case, by means of the canonical decomposition in Eq. (7), we see that 
it is possible to write any Gaussian pure state of ha + n.B modes as the direct product of ua two-mode 
squeezed vacua and ub — nA field vacua. 

Definition 4 Given a Gaussian pure state IV'g) S Q, and a bipartition into ua + ns modes, with 
nA l£ nB, the symplectic eigenvalues associated to the given bipartion are defined as 

Vk-=CQsh2rk, k = l,...,nA, (9) 
where the parameters are given by its canonical form (7). 

According to Proposition 1, a bipartition of the system induces a parameterization of the manifold of 
Gaussian states in terms of a set of non-local parameters — the symplectic eigenvalues — , and of 
a certain set of local parameters. The latter correspond to the local unitaries Uq, Uq acting on the 
subsystems, and will be denoted as a a and ub- 

All the entanglement properties of a n-mode Gaussian pure state with respect to the given bipartition 
are uniquely determined by the jia symplectic eigenvalues. The latter are invariant under local unitary 
transformations and can be explicitly computed by symplectic diagonalization [7, 18]. 

B. Main results 

The manifold Q of Gaussian states is an orbit of the group of Gaussian unitary transformations. 
Therefore, there exists a unique invariant measure on the manifold of Gaussian states induced by the 
Haar measure on the group of Gaussian unitaries. Our main interest is in computing the probability 
measure of the symplectic eigenvalues, Vk G [1, +oo), induced by that invariant measure. 

The following Theorem is the main result of this paper: 

Tlieorem 1 The invariant measure on the manifold of Gaussian pure states has the form 

UA riA 

dfiG^Kn^uA n H~^lfY['^ji''j-'^y'^~"'''^'''^t^A{aA)dtiBiaB)de, (10) 

h>k=l j=l 

where: v = {vi, . . . , Vua) '^'"^ the symplectic eigenvalues and dv = H/^ii '^^h! d^A{oiA), d^siotB) 
denote the factors in the invariant measure depending on the local degrees of freedom a a, olb, corre- 
sponding to the local Gaussian unitaries acting respectively on subsystems A, B; dd corresponds to 
the scalar term in Definition 2; and Kn.uA ^ normalization factor 
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The proof of Theorem 1 is given in Section IV A. 

This expression for the invariant measure on the manifold of ?i-mode Gaussian states makes use of a 
parameterization which highlights the symplectic eigenvalues as signatures of the entanglement across 
a given bipartition of the system. 

We also compute an explicit expression for the Haar measure on the group of Gaussian unitaries. 
For the sake of conciseness and in view of the physical applications we restrict to the subgroup of 
homogenous Gaussian unitary transformations, obtained by setting ^ = in Eq. (2). The generated 
Gaussian unitary can be written according to the Euler decomposition: 

Ug = c"'^ exp -i ^ T,ja\aj cxp ^ sual - Sk{alf exp -i ^ r/,al"j ' ^^^^ 

\ hj=l J \k=l J \ 1:3 = 1 ) 

where e M, T and T' are Hermitian matrices, and si, §2, . . . , s„ are real and non-negative pa- 
rameters (unique up to permutations). Let us recall that the Gaussian unitaries of the form Uq = 

exp iX]"j=i ^u'^i^)' ^i^h T Hermitian, define a representation of the group U(n), see e.g. [3]. 
In particular, the action of the Gaussian unitary on the canonical operators reads 



U^QakUc ^'^UkhCLh , (12) 
where U = e^^'^ £ U(n). Furthermore, 



*^-^P y X] ^J*^? ^ j exp ~ ) = cosh (2sfe) - sinh (2sfe) a[. 

Then the following Theorem holds; 



(13) 



Theorem 2 The Haar invariant measure on the group of n-mode homogeneous Gaussian unitaries 
defined in Eq. (11) is given by 



h<k=l 



where /i denotes the Haar invariant measure on the unitary group U(n), U = cxp(— iT), U' — 
exp(— iT'), A = (Ai, A2, . . . , A„), Afe = cosh2sfc, dA = 11^=1 "^^fe, and Kn is a normalization 
factor 

The proof of Theorem 2 is given in Section IV B. 



III. APPLICATIONS AND EXAMPLES 
A. A finite measure from an energy constraint 



The manifold of Gaussian states is unbounded, implying that the invariant measure cannot be nor- 
malized globally. As a consequence, all the statistical moments of the symplectic eigenvalues diverge. 
This also yields unbounded statistical moments of entropies of entanglement [2]. The same holds true 
for other entanglement quantifiers, as for instance the logarithmic negativity [27] and the coherent 
information [23], which are increasing functions of the symplectic eigenvalues. In order to avoid un- 
physical results we impose a suitable constraint, yielding an effective cutoff on the unbounded manifold 
of Gaussian states, see also [13, 18]. Here we choose to constrain the mean value of the energy in each 
subsystem. 
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Definition 5 We define the mean energies of subsystems A and B for a Gaussian pure state \'iPg)- 

^ riA 



2 



2 

A;=l 

For the sake of simplicity and motivated by the physical interpretation in the following we restrict 
to the submanifold of homogeneous Gaussian states and to the case ua = = n/2 {n even). The 
subsystems mean energies are the quantities that we will considered fixed (and finite) on a physical 
basis. The mean energies take values in the interval [1 /2, +00). They can be given a suitable expression 
by virtue of the following Lemma: 

Lemma 1 The mean energies of subsystems A andB, when calculated on homogeneous Gaussian states 
and for UA = ns ~ n/2, have the form 

^ n/2 

= 9 E \UAhk?^AhVk. (14) 



2 

h,k=l 



= \Y. \UBhk?>^Bh^k. (15) 
h,k=l 

where Ua and Ub are unitary matrices o/U(n/2), u are the symplectic eigenvalues (9), and Xa = 
{Xaii Xa2i ■ ■ ■ 1 ^An/2) ' As = (Asi, Xb2i • • • J ^Bn/2) are the parameters defined in Theorem 2. 

The proof of the Lemma is given in Section IV C. Notice that the mean energies are functions of both 
the nonlocal parameters ly and the local ones, Ua, Xa, Ub, Xb- 

By fixing the values of the local mean energies, we get an expression for the probability measure of 
the symplectic eigenvalues: 

dn{v\£A,£B) = P{v\£A,£B)<iv 

/ ^ n/2 \ / ^ n/2 \ 

= / '5 - 2 X! PAhkl'^ >^Ahl'k p - 2 E It^B/ifcP As/jJ/fc due, 

•' \ h,k=l I \ h,k=l I 



where S denotes the Dirac function. 

This integral can be evaluated by using the invariant measure of Eq. (10) after making the formal 
substitutions of d^(ayi), d^(Q;s) with d/i(ZYQ ), A^{Uq), the latter being the Haar invariant measures 
on the local groups of Gaussian unitaries. The explicit form of the Haar measures (I^IUq), d^iUQ) is 
given by Theorem 2. After integrating over all the cyclic variables, we finally obtain 



n/2 n/2 

P{v\£a,£b) = K^,n/29{y,£A)g{i^,£B) n i^^'^lfW^^^ (16) 

j<k=l 1=1 

where 

„ / ^ "/2 \ n/2 

g{iy,£A) K„/2 '^ ^^'2 ^ pAhkl'^ >^Ah'^k \ Y[ \>^Ah ~ >^Ak\dXAdii{UA) , (17) 

y h,k=l J h<k=l 

and g(z/, £b) is given by the analogous expression. 

Let us apply Eq. (16) to some simple but enlightening cases. 
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B. The case of 1 + 1 mode 

For n = 2 and ua = ns = 1, the mean energies (14)-(15) read 

£a = -^^a'^, £b — -^^Bi^ ■ 

Thus, the integral in (17) becomes 

2Ki 



Therefore, after fixing the values of the normalization constants Ki, 1^2,1 the probability measure (16) 
of the symplectic eigenvalues reads 

P{v\£a,£b) ^ —. — re ^ c 1 7' € [l,miii{£A,£B}]. 

C. The case of 2 + 2 modes 



For n — A: and ua = nB = 2, the probability measure of the symplectic eigenvalues reads 

P{vi,V2\£a.£b) = Ki,2g{vi,V2,£A) 9{vi,V2,£b) {^1 - vl)^ vlvl , 
with g given by (17). The 2x2 unitary matrix acting on the subsystem A can be parameterized as 

e''^cos(i9/2) e^^&\TL{d/2) 



Ua 



-e-"^sin(t9/2) e""^ cos (i?/2) 



with i9 G [0,7r], (p,x € [0, 27r], and the normalized Haar measure on U(2) reads dfi{UA) = 
2~i(27r)-2 sinMMtpdx- Then, from (14), (15) it follows 

f A = ^ (Aai + Xa2) {^i + ^^2) + ^ COS?? (Aai - XA2) {^i - ^2) ■ 



After integration over d<fdx we get 
K2 f 

g{vi,v2,£A) = / dAAidA/i2dcos?9|AAi - AA2I 



X(S ( - i (Aai + Aa2) {vi +1^2) -'^ cos?? (Aai - AA2) i>^i - ^2) 



2K2 



dXAl dXA2 ■ 



\Vl -'V2\ Jd 

The domain D for the A^ variables is determined by the inequality 

A£a - {Xai + Xa2){vi + V2) 



{Xai - Xa2){'>^i - V2) 



< 1 



yielding 



|l^l-I^2|[2£-(l'l+//2)]^ 



l'lV2{j^l + V2) 

Eventually, we obtain 

P{vi,v2\£a,£b) = 4i^4,2if| (i^i - V2f [2£a - {vi + y2)f [2£b - iyi + y2)f , 
where the range of the symplectic eigenvalues is determined by the inequalities 

1^1 > 1 , 

1^2 > 1, 

V2<2 min {£a, £b} ■ 

The probability measure is zero for v\ — V2, and attains its maximum when one of the two symplectic 
eigenvalues is equal to one. 
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D. Nonlinear optical parametric processes 



The measure (10) has been obtained by requiring the invariance over the whole manifold of Gaussian 
states Q. However, in physical applications, one is mostly interested in submanifolds of states arising 
from specific physical processes. The most relevant ones in which multimode Gaussian states are 
produced are nonlinear optical parametric processes, see e.g. [8]. In some of these processes there is 
a natural bipartition of the bosonic system into two disjoint subsets of canonical modes, commonly 
referred to as the signal and idler modes. Let Q (with Q C tj) be the submanifold of such Gaussian 
states. The quadratic Hamiltonians generating them are of the form 

riA "B riA nB 

Hq ^ + ^ NAhkO-AicLAk + ^ NshkasiaBk + X! X] {^'^hkaAhask + M^i^aAlasl 

h.k=l h,k=l h=l k=l 

(18) 

where Na and Nb are Hermitian matrices. Comparing with Eq. (2), the linear terms and some of the 
quadratic terms are dropped. The subgroup generated by these Hamiltonians is a projective represen- 
tation of SU(?7vi, ub) [28]. According to [25], the Gaussian states of the form Ii/jq) = cxp (— iiJ(5)|0) 
are hence coherent states for SU(n^, ub)- Following [28], these states can be written as 

( riA riB riA riB \ 

Fhk aAhask YYl ^hk a^laisl |0) . 
h=l k=l h=l k=l ) 

The canonical form for this class of states can be obtained from the singular value decomposition of 
the matrix F, Fhk = X]j ^hjfjVkj with U S U(n^), V G U(nB) and rj > are its singular values. 
By virtue of Eq. (12), we can find a pair of local Gaussian unitaries Uq, Uq such that Uq a ah I^q = 
X)j=i ^hj '^Aj and UBk^Q = ^ki ^Bi- Thus the canonical forms of these states read 

nA 

fc=i 

where the local Gaussian unitaries are of the form 

Uq = cxp -i Y TAhkO-A^aAk \ , (19) 

y h,k=l 

Uq = cxp -i Y TBhkO-BiaBk ) , (20) 

y h,k=l 

with Ta, Tb Hermitian matrices. 
The following Theorem holds: 

Theorem 3 The invariant measure on Q has the following expression 

nA riA 

d^iG = Kn.nA n {vh-VkfX{{vJ-lT''-'^^dvd^l{aA)Mc^B)Ae, (21) 

h<k=l j=l 



where Kn^uA ^ normalization constant and the local degrees of freedom a a, ctB, ore induced by the 

4^, U§. 



local unitary transformations U^, UP 



The proof is given in SubSection IV D. 

The expression for the mean energy of the subsystems takes a particular simple form on this sub- 
manifold, see Section IV C. For a balanced bipartition, ua = nB = n/2 (n even), 

^ n/2 

£a ^ £b = -^Y^i^k , 

k=l 
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and the corresponding probability measure of the symplectic eigenvalues at fixed mean energy is 

/ _^ n/2 \ n/2 

\ j=l J h<k=l 

with £ = £a = £b- 

IV. PROOFS OF THE THEOREMS 
A. Proof of Theorem 1 

In order to derive the invariant measure (10) we look at the set of Gaussian pure states as a real 
manifold. Such a manifold is embedded in the n-mode Hilbert space, which in turn is considered 
as a real vector space endowed with the real scalar product induced by the standard Hermitian scalar 
product. We then consider a system {Ha]a of generators of the algebra of the group of Gaussian 
unitary transformations, and evaluate the corresponding infinitesimal transformations on the canonical 
state; 

= -iAaU{UlHaUr)\Q) 

=: -idaUrl^aa) (22) 

where 

Ur = ^^"P [rkaAkO-Bk - rkaA\aB\j , 

k=l 



\^ha) = {UlH^Ur) |0), 

and r ~ (ri , r2 , . . . , r„^ ) G M"-* . Therefore, we can define a set of vector- valued one-forms 
{da|'0G„)}Q, where we have neglected the unitary factor —\Ur- By construction, the volume form 
generated by a maximal subset of linearly independent one-forms is the desired invariant measure on 
the manifold of Gaussian states. To compute the latter we then proceed as follows. By looking at the 
Hilbert space as a real vector space, we fix a suitable system {C,g}^i, orthonormal with respect to the 
real scalar product. We hence expand the one-forms in terms of these vectors, i.e.. 

Finally, the invariant measure is given in terms of the determinant of the matrix of coefficients: 

d/^G ~ I det J| Aq da . 

In order to obtain an expression in the form of the product of a measure on the symplectic eigenvalues 
and a measure on the local degrees of freedom, we consider a basis for the algebra of the group of 
Gaussian unitary transformations composed by: 

• the generator of the scalar-phase shift; 

• a basis for the generators which are linear in the canonical operators; 

• a basis for the quadratic generators acting on subsystem A; 

• a basis for the quadratic generators acting on subsystem B; 



10 



• a suitable set of linearly independent nonlocal generators which change the values of the sym- 
plectic eigenvalues. 

Before going into the details of the calculations, let us have a closer look at the action of the unitary 
transformation Ur on the canonical operators. The following identities hold: 

• for k < UA 

UlaAkUr = coshrfca^^. +sinhrfeaBj^ , (23) 
UlaBkUr = coshrfe asfe + sinhrfc a^^ , (24) 

• for fe > UA 

UlaskUr ^ ask- (25) 

The Lie algebra of the group of Gaussian unitaries is composed of Hamiltonians which are (at most 
quadratic) functions of the canonical operators. 

Ha = Ha{a,a^) , 

where we have introduced the short-hand notation 

t ^ t t t t ^ 

a' = (aAi,---,aA|,4,aBi,...,as;j^). 

Below we compute the corresponding vector-valued one-forms using the following relation 

M-^ha) = da {UlHa{a,a^)Ur) |0) = dai/„(a', a'^)|0) , 

where a' = U^aUr, a'^ = C/tatt/^ are expUcitly given by Eqs. (23)-(25). It follows that the vector- 
valued one-forms belong to the linear span of the following vectors: 



• for k < UA 



for k > UA 



• for h < HA , k < UA 



• for /c < UA 



• h < UA, k > UA 



\k,0,0) := a4|0), 
|0,fc,0) := asllO), 

|0,G,fc) :=asl|0), 

\hk, 0,0) -.^ GAlaAllO) , 
\0,hk,0) := asiaBilO) , 
\h,k,0) := aAiasilO) , 

|fc2,0,0) :-2-i/2(a^t)2|o^^ 

io,fc2,o) ■.= 2~'/-\aBim. 

\h,0,k) := aAlasllO) , 



\0,h,k) := aslaBllO) . 
In the following we compute the explicit expressions of the one-forms. 
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1. The scalar term 

By variation of the parameter 6, corresponding to a scalar phase-factor, we obtain the one-form 

d6'|0) . (26) 

2. The linear terms 

Let us now consider the linear generators, proportional to the complex vector ^. We consider varia- 
tions of the real parameters 

{Re{Uk)^^MUk)}k<nA , 
{Re(CBfe),Im(^i3fc)}fe<„^ , 

{Re(^Bfe),Im(^i3fc)}„^<A;<„B , 

where Re and Im denote the real and the imaginary part. Thus, using the identities (23)-(25), we get 
the following expressions for the corresponding one-forms: 

• for k < UA 

dReiUk) (coshrfc|fc,0,0) +sinhrfc|0,fc,0)) , (27) 

idlm(f^;^,) (coshrfc|fc,0,0) +sinhrfc|0,fc,0)) , (28) 

dRe(CBfc) (sinhrfe|A:,0,0) + coshrfc|0,fc,0)) , (29) 

idlm(CBfe) (sinhrfe|A:,0,0) +coshrfc|0,fc,0)) , (30) 

• for k > UA 

dRe(eBfc)|0,0,fc), (31) 

idIm(eBfc)|0,0,fc). (32) 

3. The quadratic terms: subsystem A 
We now consider the generators which are quadratic in the canonical operators of subsystem A: 

Ha^ ^ MAhk aAhaAk + MA*hk aA\aA\ + NAhk aAla^fc , 

h,k=l 

where Ma is a complex-valued symmetric matrix, and Na is a Hermitian matrix. 

Considering the variations of the parameters Re(M^^j.), lm(M^^j,), for h < k < ua, we have the 
one-forms 

dRe(M^/jfe) (cosh Tft cosh rfc|/iA;, 0, 0) + sinhrft sinhr^jO, /ifc, 0)) , (33) 
dRe(AfAfcfc)V2(coshVfc|fc2,0,0)+sinh2rfc|0,/c2,0)) , (34) 
— idIm(Myi^s,) (coshr/i coshrfe|/i/c, 0, 0) — sinhr?i sinhr^lO, hk, 0)) , (35) 
-idIm(MAfcfc)V2 (cosh Vfc|fc2, 0,0) -sinhVfe|0,fc2,0)) . (36) 

The variations of the parameters Y{.e{NAhk)^ for h < k < ua, and NAkk for k < ua, 
yield 

dRe{NAhk) (cosh r/i sinhrfc I /i, fc, 0) + sinhr/i cosh rfc I fc, /i, 0)) , (37) 

dNAkk (coshr/; sinhrfe|/c, fc, 0) + sinh^ rfcjO)) , (38) 

idl-m{NAhk) (coshrft sinhrfc|/i, k, 0) — sinhr/j coshrfe|fc, h, 0)) . (39) 
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4. The quadratic terms: subsystem B 

Moving to the local transformations on the subsystem B, we first consider the quadratic Hamiltonians 
involving the operators {ask^ 0iBk\k<nA^ which are of the form 

Hb = ^ Mshk O-BhO-Bk + Mshk O^B^asX + ^Bhk aB^aBk ■ 
hM=l 

Proceeding as in the case of subsystem A, we obtain the following one-forms from the variations of the 
matrix Mb '■ 

dRc{MBhk) (sinhr'/i sinhy-fcl/i/c, 0, 0) + cosh 7'/i cosh r/j|0, /ifc, 0)) , (40) 
dRc{MBkk) V2 {sinh^ rk\k2, 0,0) + cosh" rk\0,k2,0}) , (41) 
idlm.{MBhk) (sinhr/i sinhrfc|/ifc, 0, 0) — cosh r/i cosh rfc|0, 0)) , (42) 
idlm(MBfefe) V2(sinhVfe|A:2,0,0) -coshVfc|0,fc2,0)) . (43) 

From the variations of the matrix Nb we obtain the one-forms 

dRe{NBhk) (sinhr/i coshr/j|/i, fc, 0) + coshr/i sinhrfej/c, /i, 0)) , (44) 
dNBkk (sinhrfc coshrfcjfc, fc, 0) 4- sinh^ rfc|0)) , (45) 
—idIm{NBhk) (sinh r/i cosh rfe I ft,, fc, 0) — cosh r/i sinh rfc I A:, ft,, 0)) . (46) 

Then we consider the Hamiltonians coupling the operators {ask: o,Bk}k<nA with {asj,, aBj^}fc>„^, 
that is, 

riA riB 

= ^ ^ {PB'hkaBhaBk + PB'hk^'siaBl 

h—1 k—nA-\-^ 

+QB'hk (^Bh^B^k + Q*B'hk ^B^aBl^j , 

where Pb' and Qb' are complex-valued matrices. For h < ua and ua < k < ub get 

dRc(PB,,,fe) coshr,,|0,ft,fc), (47) 

idIm(PB,,,fc) coshr,,|0,ft,fc), (48) 

dReiQB'hk) sinhr,,|ft,0,fc), (49) 

idlm(QB'ftfc) sinhrft|ft,0, k) . (50) 

Finally, let us consider the quadratic Hamiltonians containing only the operators {as^, as[.}fc>„^: 

rtB 

Hb' = X] ^'^B'hk aBhO-Bk + MB'lk aslasl + ^'B'hk aslciBk , 

h.k—nA-\-^ 

where the matrix Mb' is complex-valued and symmetric and the matrix Nb> is Hermitian. We obtain 
that the variations of the elements of the matrix Nb' yield vanishing one-forms. The non-zero one- 
forms are generated by the variations of Mb', 

dReiMB'hk)mO,hk), (51) 

dRe(Ms'fcfc)y2|0,0,fc2>, (52) 
-idIm(MB,;,.fe)2|0,0,ftfc), (53) 

-idIm(il/B,fcfc)V2|0,0,fc2). (54) 
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J. The nonlocal generators 

It remains to consider the nonlocal transformations whose action changes the symplectic eigenvalues. 
These are generated by Hamiltonians of the form 

{HNL)k = i {aAkask ~ aAlasl) , k < ua , 
The corresponding one-forms are 

-idrfe|fc,fc,0) = -^^^^|fc,fc,0), k<nA. (55) 
2 smh 2rfe 

6. The invariant measure 

We can now compute the invariant measure on the manifold of Gaussian states. 

First, we consider the one-forms corresponding to linear Hamiltonians. The matrix of coefficients 
is readily obtained from Eqs. (27)-(32), from which we obtain the following factor in the invariant 
measure 

ti=l k=l 

Second, we consider the one-forms corresponding to the scalar phase-shift and the quadratic Hamil- 
tonians. The matrix of coefficients can be straightforward obtained from Eqs. (26), (33)-(39), (40)-(46), 
(47)-(50), (51)-(54), and (55). A maximal subset of linearly independent one-forms can be obtained by 
eliminating the one-form in (38), which is proportional to the one in (45) due to the symmetry of the 
canonical form (7). 

From the matrix of coefficient one gets 

TIA riA 

idetj|=c n K--.^)'n-i(-i-i)^"™'' 

h<k=l j=l 

where C is a constant factor. Finally, the invariant measure in Eq. (10) is obtained by inserting the 
differentials d9, di/i , . . . , dvnA ' identifying the factors depending on the local degrees of freedom 



dAi(aA) = dU n d^MAhkd^NAhklid^MAu, 

h<k=l i=l 
riB riA 

d^i{aB) = die n dHIBhkd^NBhkX{dHIB^^dNB, 



■II 

h<k=l i=l 
iiA riB riB 

n ^^PB',id?PB'jid^QB',id^QB',i n ^^^^b-^- 

j=ll=nA + l p<g=nA+l 



B. Proof of Theorem! 

In order to derive the explicit expression of the Haar measure on the group of n-mode homoge- 
neous Gaussian unitaries, we apply the unitary transformations on the vacuum state. Using the Euler 
decomposition in Eq. (1 1) we obtain the n-mode homogeneous Gaussian state 

IV'g) = Ug\0) = e-'' exp I -i T^,aU, j cxp 8^4 - Sk{alf^ |0) , 



14 



where we have used cxp ^— i ^/j^l^j) 1^^) = l^^)- Then, by variation of the parameter 9, of the 

elements of T and of the parameters s ~ (si, S2, ■ • • Sn), we obtain a set of vector valued one-forms. 
By proceeding as in Section IV A, these one-forms can be used to derive an explicit expression for the 
invariant measure on the manifold of homogeneous Gaussian states. Then, the Haar measure on the 
group of homogeneous Gaussian unitaries can be readily obtained from the latter. 
By variation of the parameter 9, we obtain the vector-valued one-form 

d6i|0). (56) 

The variations of the parameters s yield 

idskV2\k2), (57) 
where |/c2) = 2^^/^(a[)^ |0). The variations of the parameters Tkk yield the one-forms 



kk 



(58) 



2-i/2sinh2sfe|A:2) + (sinhsfc)2|0) . 

Similarly, by variations of the parameters 'Re{Thk), we obtain 

dRe{Thk) sinh (s,, + Sh)\hk) , h<k<n, (59) 
and the variations of the parameters Im(Tfik), yield 

idlm{Thk) sinh (sfe - Sh)\hk) , h < k <n, (60) 

where \hk) = a\al\0). 

The volume form generated by these one-forms is by construction the invariant measure on the 
considered submanifold of Gaussian states. First, we notice that the volume forms generated by the 
variations of the matrix elements of T and T' are 

n n 

Y[dTkk n d'^^J = dAi(C/), 

n n 

HdT^,. n = dM(C/'), 

k=l i<j=l 

where d/i(C/) and d/i(C/'), with U ~ cxp(— iT) and U' = cxp(— iT'), denote the Haar measure on the 
unitary group \J{n). Then, from the one-forms (56), (57), (58), (59), (60), we are able to derive the 
following expression for the Haar measure on the group of n-mode homogeneous Gaussian unitaries: 



dM(Z^G)- n \^h-Xk\Y[d\,dfi{U)dfi{U')^ 

h<k=l 3=1 



with Xk = cosh 2sk- 



C. Proof of Lemma 1 

For a given rt-mode homogeneous Gaussian state |i/'g)i and a bipartition of the system defined by 
two disjoint sets of canonical operators 

{aAk,CiAl}k=l,...nAJ {asfc, asl.}fe=l,...ni5 , 

riA + TT-B = n, we consider the mean value of the energy of one of the two subsystems. The non- 
homogeneous case can be analyzed in a similar way, and gives rise to an additional term in the mean 
energy. To fix the ideas we consider the mean energy of subsystem A, 

^ riA 

= ^ ^{^g\ (aAittAk + aAkO-Ai^j IV'g) , 
fc=l 
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and restrict to the case of balanced bipartition, ua = ns = n/2. 
Using the canonical form (7) we get 

^ n/2 

" k=l 
^ n/2 

" k=l 

where a^^. = Uq^ aAkl^c '^aI ~ ^0^*^-41^0 linear combinations of the operators 
{aAk, aAl.}A:=i,...m- Their explicit form can be written starting from the Euler decomposition of the 
homogeneous Gaussian unitary, see Eq. (11), 

/ n/2 \ /n/2 \ / n/2 

Uq = e^^'^cxp -i ^ TA^JaAlaAJ exp ^SAkO-Al - SAk{aAl)^ exp -i ^ TA^jaAlaAj 

\ = l / \fc=l / \ tj=l 

Then, Eqs. (12), (13) imply 

^ n/2 

'S'a = 2 5Z l^Ahkl"^ ^Ah l^k , 
h,k=l 

where Ua = exp (— ITa), Xau = cosh2sAk., and z^^'s are the symplectic eigenvalues. The analogous 
expression is obtained for the mean energy of subsystem B, 

^ n/2 

'S's = 2 X/ Wshk?' ^Bhl^k ■ 
h,k=l 

The local mean energies take a particular simple form for the submanifold of states considered in 
Section HID. In that case, from Eqs. (19), (20) we get Xau = ^Bk = 1, which in turn implies 

^ n/2 

£a ^ £b = -j^ ^ Vk ■ 

k=l 

D. Proof of Theorem 3 

To derive the invariant measure in Eq. (21), we proceed along the same steps of Section IV A, with 
the difference that only the terms which are compatible with the form of the Hamiltonian in Eq. (18), 
and hence compatible with the local Gaussian unitaries in Eqs. (19), (20), have to be retained. Thus, 
the invariant measure is generated by the one-forms (26), (37)-(39), (44)-(46), (49)-(50), and (55). 
This yields the expression in Eq. (21) where the factors depending on the local degrees of freedom are 
explicitly given by 

n-A 

d/x(aA) = Jl d^NAhk, 

h<k=l 

riB riA riA "b 

dK^B) = n 'i^NBhkl[dNBull n d^QB'.id^QB'ji- 
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